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Abstract 

The inclusive production of prompt photons with energy above 10 GeV is measured using 
the OPAL detector in hadronic Z° decays at LEP. In contrast to previous measurements, the 
prompt photons were not required to be isolated. The production rate and energy spectrum 
are found to be in agreement with QCD predictions for the quark-to-photon fragmentation 
function. 
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1 Introduction 



We gain a greater understanding of the properties of the elementary building blocks of matter 
and their interactions by studying the properties of hadrons, leptons and photons produced 
as a result of a primary interaction. The properties of hadrons and charged leptons produced 
in e + e~ collisions have been studied in great detail at different centre-of-mass energies. In the 
case of photons radiated off quarks, prompt photons, in hadronic e + e~ collisions, much less 
information is available due to the difficulty of separating these photons from those produced 
in the decays of other particles [|IJ]-|§. Both the shape and normalisation of the inclusive 
prompt photon energy spectrum in e + e~ collisions are predicted, through the calculation of 
the quark-to-photon fragmentation function, by leading-order perturbative QCD || [7|]. This 
asymptotic prediction has been parametrised in ||. Non-perturbative effects can be included 
in the calculation through the vector-meson dominance ansatz as in |], ID| , where boundary 



terms missing in || were also accounted for. The higher-order terms were calculated, and 
seen to be small at the energies close to the Z° peak. Direct experimental study of these 
predictions is important in providing insight into the non-perturbative and higher order effects 
in the radiation of photons from quarks. This will also make theoretical predictions of photon 
production in other processes, such as those occuring at pp and pp colliders, more reliable, 
thus improving sensitivity to possible new phenomena. 

At LEP the first measurement of prompt photon production in hadronic Z° decays was 
made by the OPAL Collaboration [JTTJ] for photons isolated from other particles in the event, 
as suggested in []T2[. The production of isolated prompt photons was studied in great detail 
by all LEP experiments [fL3[j-||16|. Following the suggestion of [0, the ALEPH Collaboration 
extracted the quark-to-photon fragmentation function from the study of non-isolated photons 
in jets containing a photon carrying more than 70% of the jet energy [|1 



Here we present a measurement of the inclusive prompt photon energy spectrum in hadronic 
Z° decays at LEP. This method of studying the quark-to-photon fragmentation function was 
suggested in To separate prompt photons from the photons from decays of other 

particles we use the following method. We selected clusters in the electromagnetic calorimeter 
not associated with charged tracks. A set of cuts were applied to reduce the background in 
the sample. The distribution of a variable characterising the transverse shape of the clusters 
in data was then fitted with a linear combination of the distributions for photons and for 
background to determine the fraction of prompt photons in the selected sample. The result 
was then corrected for the selection efficiencies, detector effects and initial state radiation. In 
the following sections we describe the OPAL detector, the event and electromagnetic clus- 
ter selection (sec. and the determination of the number of photons in the selected sample 
(sec. H). The efficiency and acceptance corrections are described (sec. |5]) followed by the 
study of systematic effects (sec. ^]). Finally the measured prompt photon energy spectrum is 
presented and discussed (sec. [F]). 
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2 The OPAL detector 



The OPAL detector operates at the LEP e + e~ collider at CERN. A detailed description of the 
detector can be found in pPfl . For this study, the most important components of OPAL were 
the central detector and the barrel electromagnetic calorimeter with its presampling detector. 
The central detector, measuring the momenta of charged particles, consists of a system of 
cylindrical tracking chambers surrounded by a solenoidal coil which produces a uniform axial 
magnetic field of 0.435 T along the beam axis]]. The detection efficiency for charged particles 
is almost 100% within the polar angle range | cos^l < 0.95. 

The electromagnetic calorimeters completely cover the azimuthal range for polar angles 
satisfying | cos#| < 0.98 providing excellent hermeticity. The barrel electromagnetic calorime- 
ter covers the polar angle range | cos#| < 0.82. It consists of 9440 lead glass blocks, each 24.6 
radiation lengths deep, almost pointing towards the interaction region. Each block subtends 
an angular region of approximately 40 x 40 mrad 2 . Half of the block width corresponds to 
1.9 Moliere radii. Deposits of energy in adjacent blocks are grouped together to form clus- 
ters of electromagnetic energy. The intrinsic energy resolution of a E /E = 0.2% © 6.3%/ VE 
is substantially degraded (by a factor ~ 2) due to the presence of two radiation lengths of 
material in front of the lead glass. For the intermediate region, 0.72 < | cos#| < 0.82, the 
amount of material increases up to eight radiation lengths causing further degradation in the 
energy resolution. The two endcap calorimeters, each made of 1132 lead glass blocks, 22 
radiation lengths deep, cover the region of 0.81 < | cos#| < 0.98. In this study the measure- 
ment of inclusive photon production is restricted to the barrel part of the detector. Most of 
the electromagnetic showers start before the calorimeter and their position at the entrance of 
the calorimeter is measured by a barrel electromagnetic presampler made of limited streamer 
mode chambers. The presampler covers the polar angle range | cos#| < 0.81 and its angular 
resolution for photons is approximately 2 mrad. 



3 The selection of events and electromagnetic clusters 

Our study was based on a sample of 2.5 million hadronic Z° decays selected as described in 



|2l| from the data accumulated with the OPAL detector at LEP in 1992, 1993 and 1994 at an 
e + e~ centre-of-mass energy of 91.2 GeV. We did not use the off-peak data to avoid additional 
complications in dealing with data collected at different e + e~ centre-of-mass energies. We 
required that the central detector and the calorimeters were fully operational. Temporary, 
local inefficiencies in the presampler chambers were monitored and taken into account. 

To study the properties of the background we used Monte Carlo events produced with 
the parton shower generators JETSET 7.4 @ (3.9 million events) and HERWIG 5.8 [0(1.1 



million events) with generator parameters given in [24]]. We also used samples of events with 



only single photon or a 7r° meson present in the detector. The Monte Carlo (MC) samples 

1 In the OPAL coordinate system the x axis points towards the centre of the LEP ring, the y axis points 
upwards and the z axis points in the direction of the electron beam. The polar angle 6 and the azimuthal 
angle 4> ar e defined with respect to the z and x-axes, respectively, while r is the distance from the z-axis. 
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were passed through the full simulation of the OPAL detector |25| and subjected to the same 
reconstruction and analysis procedure as the data. 

The difficulty in the measurement of prompt photon production lay in the separation 
of the signal from background. The QCD shower models predicted a signal-to-background 
ratio of approximately 1/200. Background clusters, with no charged track associated with 
them, are dominated by photons from decays of hadrons, particularly 7r° — > 77 (~ 57%) and 
7] — > 77 (~ 10%). The other sources of background, like the interaction of neutral hadrons 
such as Kl's or neutrons in the material of the calorimeter contribute at the level of few per 
cent each. More than one particle can also contribute to a cluster. The transverse profile 
of the calorimeter cluster can be used to differentiate between clusters coming from different 
sources. The hardest to remove are clusters produced by it — > 77 and rj — > 77 meson decays 
because they can be very similar, especially for higher cluster energies, to those produced by 
one photon. An irreducible, but very well predicted, background comes from the initial state 
radiation (ISR) photons radiated by the beam particles before they interacted. 



The tracks and calorimeter clusters were selected as described in |21]]. In addition we 
required a cluster energy to be larger than 10 GeV and cluster polar angle such that | cos#| < 
0.72. We then applied three cuts motivated by studies with simulated events. 



Cut 1. There was required to be no charged track associated with the cluster. Tracks were 
extrapolated to the calorimeter surface. A track was associated with a calorimeter 
cluster if it extrapolates to the calorimeter within 24 mrad (approximately half of a 
lead glass block width) of the centre of gravity of a cluster. In Figure [TJa and [IJb, the 
normalised distributions of the angle between the calorimeter cluster and the nearest 
track are shown for background clusters and prompt photons in the JETSET model, 
for small and large cluster energies. For the higher energy clusters a contribution from 
electrons and positrons from conversions of prompt photons in the beam pipe and central 
detector is seen for small angles 5. For lower energy clusters this effect is diluted both by 
the greater separation between the extrapolation of the track to the calorimeter surface 
and the centre of gravity of the calorimeter cluster, and by the presence of other tracks 
in the proximity of the prompt photon. 

Cut 2. We required the presence of a presampler cluster within 24 mrad of the centre of gravity 
of the calorimeter cluster. The differences in azimuthal \A<f>\ and polar |AG| angles 
between the positions of the calorimeter cluster and the presampler cluster were required 
to satisfy A = min(|A0|, |A0|) < A cut . The distributions of A for single, isolated 
photons and 7r°'s in the detector as well as all background clusters in JETSET events 
are shown in Fig. 0a and |2|b. The value of A tended to be larger for clusters produced 
by two overlapping photons, e.g. from a 7r° — ► 77 than for clusters produced by a 
single photon, because the presampler measures the cluster position at an early stage of 
development of the electromagnetic shower. 

Cut 3. The transverse profile of the cluster in the calorimeter was required to be compatible with 
that produced by an isolated photon. It was first assumed that the cluster was produced 
by a photon. The impact point of the photon was varied until the best description of the 
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observed lateral shower profile by a reference profile was found. The reference profile was 
obtained by the parametrisation of the results of MC simulation of the isolated photon in 
the detector. The fast algorithm described in section 4.1 of [26] was used. The resulting 
variable, S, is proportional to the \ 2 f° r matching the measured and predicted energy 
sharing between the calorimeter blocks. The distributions of S for single photons, single 
7r°'s and background clusters in JETSET events are shown in Fig. [|c and Qd. 



The energy dependent values of A cut and S cut were chosen to optimise the separation of signal 
from background and are shown in Tab. |]. 

In total 23106 clusters passed the selection procedure. The signal-to-background ratio, 
estimated from simulation, was improved from 1/200 to 1/130 after cut 1, then to 1/50 after 
cut 2 and finally to 1/6 after cut 3 with respectively 92%, 51% and 42% of the signal retained. 
According to the simulation the background clusters passing the selection were produced 
mostly by 7r° — > 77 (~79%) and rj — > 77 (~11%) decays. 

The cuts lead to a strong reduction of the background while prompt photons were much less 
affected. The fraction of photons rejected by the cuts can be corrected for and its knowledge 
will affect the systematic uncertainty of the measurement as detailed in Sections || and |[ 
The efficiency of cuts 2 and 3 for photon clusters well separated from other particles in the 
event was determined directly from the data using a sample of photons in radiative lepton 
pair events e + e~— > £ + £~7 (£ — e, //). In addition, in hadronic events, further losses of prompt 
photons occur when other particles hit the calorimeter close to the photon. This can result in a 
cluster being associated to a track or being sufficiently distorted to fail the selection criteria. A 
small fraction (~ 6%) of prompt photons also converted in the beam pipe or central detector. 
The correction for these effects was estimated using Monte Carlo as detailed in sec. [S[ In 
Figure [I] we compare the data and MC distributions for the angle 5 between the cluster and 
track closest to it. Differences between data and Monte Carlo are concentrated in the region 
of small angles, below our cut value of 24 mrad, a region containing only a small fraction of 
prompt photons. In Figure ||| we present data and MC distributions for the variables S and 
A used in cuts 2 and 3. We compare distributions for photons from radiative lepton events 
with MC for single, isolated photons (a and c). We also show distributions for data and MC 
clusters from the t + t~ events (r ± — > p^v T and p ± — > 7r ± 7r , section 4.2 of (b and d). The 
differences between data and MC are small and concentrated mostly in the region below our 
cut values, so they will cause only small systematic effects. 



4 The determination of the number of photons in the 
selected sample 

We determined the fraction of photons in the sample remaining after cuts 1-3 above using the 



cluster shape fit variable C used in the previous OPAL studies of photon production |TT], |16 
The fit algorithm applied was more sophisticated than that used in the the calculation of 
the variable S. The variable C had a better background rejection power than S. Due to 
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the similar fit algorithms and shower parametrisations the C and S variables are correlated, 
although not fully. The definition of C is 

N b Y K red ) 2 + « bs ) 2 ' 

where: _E° bs is the energy observed in calorimeter block number i; i?f red is the predicted energy 
in calorimeter block number i; af md and <r° bs are the energy dependent errors on Ef red and 
E° bs , respectively; and iV 6 is the number of blocks in the cluster. Ef md was taken from the 
best fit of the shower profile parametrisation, assuming that the cluster was produced by a 
isolated photon, to the observed energy sharing between the calorimeter blocks. The reference 
profiles varied as a function of cos 9 because of the varying amount of material in front of the 
calorimeter. 

We fitted the distribution of C in the data with a linear combination of MC distributions 
for photons and background for clusters passing the same selection criteria as the data: 

C fi t = fC, + (1 - f)C hkg . (2) 



The fraction / of photons in the selected sample was the fit parameter. For the distribution 
of the background, Cbkg, clusters from JETSET hadronic Z° events were used, where initial 
state radiation and prompt photons were removed from the sample. For the distribution of the 
photons, C 7 , we used a simulated sample of isolated photons. The fit for / to the C variable 
distribution in the data was performed separately in seven bins of cluster energy as shown in 
first column of Tables [l| and |[ The C 7 and Ck g distributions had only a small dependence 
on cluster energy within a given bin. A binned maximum likelihood method [27 was used to 
fit the C variable distribution between and 5. 

Since we used simulated distributions for C 7 and Cbkg in eq. |2|, it is crucial to check that 
the C variable is well described in the simulation. We show data and MC distributions for 
photons (from radiative lepton pair events and single, isolated photon MC) in Fig. |3|e and 
clusters from data and MC r decays (r ± — ► p ± z/ r and p ± — > 7r ± 7r , sec. 4.2 in |26|) in Fig. |3|f. 



The simulation describes well the C variable distributions for clusters produced by isolated 
photons as well as by 7r°'s. In Fig. |5| we show distributions of the C variable for background 
clusters from different sources in the JETSET simulation. The shapes of the distributions are 
similar, although tt° — > 77 and rj — > 77 decays tend to produce less clusters with higher values 
of C than other sources of background. 

The fitted fraction of photons / is shown in Tab. |2| for different cluster energy ranges. 
The comparison of the data and fit results is given in Fig. ^. The contributions from prompt 
photons and background to the fit are shown. The x 2 \ taking into account the statistical 
errors on the C 1 and Cbk g distributions, were between 16 and 37 for 23 degrees-of-freedom. 
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5 Corrections for efficiency, acceptance and initial state 
radiation 



The energy spectrum of prompt photons obtained in the previous section was corrected for 
photons lost in the selection process and outside the geometrical acceptance. To ensure 
that the energy spectrum of photons was not biased, efficiency corrections were determined 
separately for each energy bin. Then the contribution due to initial state radiation was 
subtracted and the corrected photon energy spectrum was normalised to the total number of 
hadronic events. 

We applied corrections for the following effects. 

1. Local, temporary inefficiencies in the presampler system. This factor was determined 
from data to be 1.23 with negligible statistical error. 

2. Rejection of real photons by cuts 2 and 3. A correction was determined using the data 
sample of photons in radiative lepton pair events. 

3. Additional rejection, by the combined effects of cuts 1, 2 and 3, of photons that formed 
calorimeter clusters with other particles in the event or converted in the beam-pipe or 
central detector. This correction, called photon environment in Tab. [1], was determined 
with the JETSET Monte Carlo as the ratio of the combined cut efficiency for prompt 
photons where no other particle contributed to the calorimeter cluster, to the efficiency 
for all prompt photons. 

4. The contribution of initial state radiation, estimated using the KORALZ program ver- 
sion 4.0, was subtracted. 

Values of the energy dependent corrections are shown in Tab. [I]. To compare our result with 
theoretical predictions we applied an additional correction for the rejection of photons by the 
cut on the polar angle 9 {\co§6\ < 0.72). The correction was 1.58. We assumed the leading- 
order 1 + cos 2 9 dependence of the photon production cross-section, which was consistent with 
the polar angle distribution of the selected calorimeter clusters. The fully corrected energy 
spectrum of prompt photons in hadronic Z° decays is shown in Tab. ^| and in Fig. || 

6 Systematic effects 

We checked the dependence of our result on possible deficiencies in the simulation of the de- 
tector and on the particle composition of the background. The main sources of the systematic 
uncertainties were estimated as follows: 

1. The sensitivity of the fit result to the quality of the MC reproduction of the C variable 
was determined as follows, separately for each cluster energy bin. The C 1 and Cbk g 
distributions used in the fit were simultaneously scaled by (1 ± a), where a (~ 4%) was 



9 



the error on the mean value of C from the Monte Carlo added in quadrature to that 
from the data. These modified distributions were then used in the fit. The differences 
between these results and those obtained with unmodified distributions were assigned 
as the systematic error. These uncertainties, ranging from 10-23% and partly correlated 
between different photon energy points, were the dominant source of systematic error in 
our measurement. 

2. The fractions of photons in the data obtained from the fit using the background spectrum 
predicted by the HERWIG MC were in agreement with those obtained using the JETSET 
MC background. This showed that, within our statistical precision, the result did not 
depend on the details of the model implementation of the parton shower development 
and hadronisation processes in the MC generators. The difference between the fraction 
of prompt photons obtained from the data fitted with the Cbkg predicted by the JETSET 
and HERWIG was assigned clS 8u systematic error. 

3. Although we relied on the data as much as possible in the determination of the efficiency 
corrections we had to resort to MC to estimate how the efficiencies for photons passing 
the selection cuts were modified by the presence of other particles in the event (correction 
3 in sec. |5|). The corrections obtained with JETSET and HERWIG models gave results 
consistent within statistical errors despite the differences in the modelling of the prompt 
photon radiation, parton shower and hadronisation for the two models. The difference 
between corrections obtained with the JETSET and HERWIG models was assigned as 
a systematic error. 

In addition several further checks were performed, none of which produced a statistically sig- 
nificant, at the one standard deviation level, difference from the result of the default procedure 
and were not included in the systematic error. 

• In equation ||] the amount of the background in the sample was the fit parameter. There- 
fore, the fitted value of / was not, to first order, sensitive to the background flux as 
incorporated in the MC generator. In principle, some sensitivity to the relative fluxes of 
different background sources remained, since it could change the shape of the C variable 
distribution of the background. This sensitivity was estimated by repeating the fit with 
the number of clusters in the background produced by rj mesons, the second largest 



source of background, adjusted by the uncertainty on the rj yield of ±20% [^9j. This 
resulted in 1.5% change of the prompt photon yield. 

The systematic error from the modelling of the C variable (point 1) was consistent with 
an estimate which assumed a dependence of the factor a on the value of C variable. For 
C greater than 2 the factor a was put to zero. For lower C values a linear dependence 
was assumed, such that the mean value of a for C between and 2 was equal to the 
value of a used in the default procedure. 

The fit of equation |2| was repeated with the distribution of C for prompt photons in 
JETSET used as C 7 instead of the default of isolated photons. 

The values of the cuts A cut , S cut and for the association cut were changed by ±20 %. 
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The fit result (section |]) did not depend on the bin size of the fitted distribution (changed 
by a factor of 2) or the upper fit boundary (moved between 5 and 10 in C). 

The effect of the energy resolution was negligible for the size of the photon energy 
bins used. The absolute energy scale for the electromagnetic calorimeters is very well 
calibrated using Bhabha scattering events and cross-checked with, for example, the ir° 
mass peak position. 

The background of calorimeter clusters from r + r~ events in the selected sample (esti- 



mated with KORALZ 4.0 [23) was less than 0.5% of the observed photon signal. 



The contamination from clusters produced by the LEP accelerator background or cosmic 
rays, estimated with the events collected using a random beam-crossing trigger, was 
below the level of 2 clusters per million events and no cluster passed the selection criteria. 



To summarise, the principal systematic error comes from the uncertainty on the fitted 
fraction of photons in the data due to the uncertainty on the quality of the MC reproduction 
of the C variable (point 1). The uncertainties due to the possible dependence of our result on 
the modelling of the prompt photon radiation, parton shower and hadronisation in the MC 
generators were estimated in points 2 and 3. The total systematic uncertainty was calculated 
as the sum in quadrature of the propagated errors listed in points 1, 2 and 3, as shown in 
Tab. |||. The systematic errors were partly correlated between the different photon energy 
ranges. 



7 Summary and discussion 

The energy spectrum of prompt photons in hadronic Z° decays, corrected for the geometrical 
acceptance as discussed in section [|, is shown in Tab. |2| and in Fig. |6|. The inner error bars 
are statistical, and the outer combined statistical and systematic (added in quadrature). The 
breakdown of the systematic uncertainties is given in Tab. || 

To compare our result with leading-order theoretical predictions for the quark-to-photon 
fragmentation function we use the leading-order cross-section for prompt photon production 
in e + e~ annihilation given by formula (12) in |Tj| : 

1 da(E 7 ) 4 _ 

w i D y/<v ( 3 ) 



Chad dE 1 y/s 



q 



where E 1 is the photon energy; is the e + e~ centre of mass energy; w q is the relative con- 
tribution from quark flavour q (w q = r z _> q q/T z _+ hadrons ); o"h a d is the production cross-section 
for e + e~ hadronic events, and -D 7 / q is the quark q to photon fragmentation function. In Fig. |6| 
we plot the QCD prediction using an asymptotic leading-order D y / q |6|, [7J as parametrised 
by Duke and Owens g with Q 2 = Mf and A = 0.2 GeV. In Fig. | we also plot the QCD 
predictions using D y / q from leading-order (LO) calculations by Gliick, Reya and Vogt ||. We 
also show the higher-order (HO) prediction for the prompt photon production cross-section 
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by Gliick, Reya and Vogt M and the beyond leading logarithm (BLL) prediction by Bourhis, 
Fontannaz and Guillet flTOfl . The factorisation schemes used were DIS 7 and MS respectively. 



The authors of M and JTDJ included in their calculations non-perturbative effects through the 



vector-meson dominance ansatz, although using different experimental inputs. 

Our data are in agreement with these theoretical predictions. The experimental precision 
is not sufficient to discriminate between them. The ALEPH data on the production of jets 
containing a photon carrying a substantial fraction (above 70%) of the jet energy (Fig. 4 in 
|T8| ) show clear disagreement with the Duke-Owens parametrisation M. It was noted in ]10 



that a possible reason is that the ALEPH measurement is restricted to only part of the phase- 
space for photon production in hadronic Z° decays. 

To summarise, we have measured the inclusive production of prompt photons with energy 
above 10 GeV in hadronic Z° decays. Good agreement is found with current QCD predictions 
for the quark-to-photon fragmentation function. 
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Table 1: Energy dependent cuts and corrections to the photon energy spectrum. Corrections 
for photon efficiency and environment are multiplicative, ISR is subtracted. 



E, GeV 


mrad 


<Scut 


isolated photon 
efficiency 


photon 
environment 


iV^/GeV/lO 6 
events 


10-15 


3. 


1.5 


2.51 ±0.13 


1.34 ±0.10 


6.71 ±0.13 


15-20 


4. 


2.5 


2.06 ±0.12 


1.36 ±0.11 


2.62 ±0.08 


20-25 


5. 


4. 


1.51 ±0.07 


1.37 ±0.10 


1.52 ±0.06 


25-30 


5. 


5. 


1.47 ±0.08 


1.25 ±0.09 


1.27 ±0.06 


30-35 


5. 


5. 


1.42 ±0.08 


1.20 ±0.09 


1.28 ±0.06 


35-40 


5. 


5. 


1.60 ±0.09 


1.14 ±0.10 


1.20 ±0.06 


40 - 45.6 


5. 


5. 


1.57 ±0.05 


1.03 ±0.13 


2.87 ±0.09 



Table 2: The fraction / of prompt photons in the selected sample and the corrected energy 
spectrum of prompt photons in hadronic Z° events (the first error is statistical, the second 
systematic). 



E, GeV 


Mean E, GeV 


/,% 


iV 7 /GeV/10 6 events 


10-15 


12.0 


32.7 ± 2.1 


625 


±70 


±99 


15-20 


17.3 


18.2 ± 1.7 


266 


±36 


±63 


20-25 


22.3 


12.9 ± 1.4 


142 


±20 


±37 


25-30 


27.3 


27.9 ± 1.9 


212 


±24 


±39 


30-35 


32.3 


25.7 ± 2.8 


118 


±17 


±28 


35-40 


37.1 


29.9 ± 5.0 


75 


±15 


±20 


40 - 45.6 


42.1 


80.4 ± 10.2 


48 


± 9 


± 8 
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Table 3: Contributions to the systematic error on the photon energy spectrum shown in Tab. 2. 



Mean E, GeV 


iV 7 /GeV/10 6 events 


Systematic errors on iV 7 /GeV/10 6 events 


description of 
C variable 


model dependence 
of background 
distribution 


photon acceptance 


12.0 


625 ± 70 ± 99 


±60 


±73 


±29 


17.3 


266 ± 36 ± 63 


±54 


±31 


±12 


22.3 


142 ± 20 ± 37 


±32 


±17 


± 6 


27.3 


212 ± 24 ± 39 


±28 


±25 


±10 


32.3 


118 ± 17 ±28 


±23 


±14 


± 6 


37.1 


75 ± 15 ± 20 


±17 


± 9 


± 4 


42.1 


48 ± 9± 8 


± 5 


± 6 


± 3 




Figure 1: The angle, 5, between the calorimeter cluster and the nearest track for hadronic 
events: for the data, background clusters in the JETSET simulated event sample and 
prompt photons in the JETSET sample, for different cluster energies: 10 < E < 15 GeV 
a), and 35 < E < 40 GeV b). Distributions are normalised such that the total yield for 
each curve is unity. Most of the photon signal is off scale, with 5 > 100 mrad. The regions 
removed by the cut are shown by arrows. 
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Monte Carlo simulation 
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Figure 2: The distributions of the variables A (a and b), and S (c and d) for clusters 
from simulation of single isolated photons and 7T°'s and for all background clusters in 
JETSET events (after cut 1), for different cluster energies: 10 < E < 15 GeV a) and c), 
and 35 < E < 40 GeV b) and d). Distributions are normalised such that the total yield 
for each curve is unity. The regions removed by the cuts are shown by arrows. 
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Figure 3: The distributions of the variables A (a and b), S (c and d) and C (e and f) 
for photons in radiative lepton pair data events and simulated single, isolated photons 
(a, c and e); and 7r°'s in — > p ± u T , p* 1 — > 7r ± 7r° decays (b, d and f). Distributions are 
normalised such that the total yield for each curve is unity 
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Figure 4: Comparison of the data and fit for different cluster energy ranges. Contributions 
from prompt photons and background in the fit are shown. 
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Monte Carlo simulation 
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Figure 5: C variable distributions for background clusters from different sources in the 
JETSET simulation for different cluster energy ranges. 
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Figure 6: The photon energy spectrum in hadronic Z° decays compared to various theo- 
retical predictions: the Duke-Owens parametrisation [8], the Gliick, Reya and Vogt pre- 
dictions including leading-order (LO), higher-order (HO) and higher-order without the 
non-perturbative corrections (HOpl) [9]. The Bourhis, Fontannaz and Guillet prediction 
shown include effects beyond leading logarithms (BLL) [10]. 
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AXODRAW 



J . A . M . Ver maser en 

NIKHEF-H 
P.O. Box 41882 
1009 DB Amsterdam 



Abstract 

Axodraw is a set of drawing primitives for use in M^jX. These can be 
used for the drawing of Feynman diagrams, flow charts and simple graph- 
ics. Because it uses postscript for its drawing commands it works only 
in combination with the dvips of Radical Eye Software which is presently 
the most popular dvips program. More will be added in the future. It 
allows whole articles including their pictures to be contained in a single 
file, thereby making it easier to exchange the article file by e-mail. 
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1 Using Axodraw 



The file axodraw. sty is a style file for I^TeX. It should be included in the 
documentstyle statement at the beginning of the document. An example would 
be: 

\documentstyle [a4 , 1 lpt , axodraw] {article} 

Because axodraw. sty reads also the epsf.sty file that comes with many imple- 
mentations of TgX and in particular those that rely on the dvips program by 
Radical Eye Software for the printing, this file should be present in the system. 
If this file is not available one should obtain it from another system. The author 
feels in no way responsible for the problems that may occur when a different 
dvi-to-postscript program is used. 

The drawing is actually done in postscript. Because the above mentioned 
dvi-to-postscript converter allows the inclusion of postscript code the graphics 
primitives have been included in the file axodraw. sty in terms of postscript. If 
another postscript converter is used, one may have to adapt the syntax of the 
inclusion of this code to the local system. 

The commands of Axodraw should be executed inside either the picture or 
the figure environment. Inside this environment it is possible to place objects 
at arbitrary positions and put text between them. In principle one could try 
to draw objects with the facilities of WTpJi itself, but it turns out that the 
commands in the picture environment are not very powerful. Axodraw gives 
good extensions of them. An example would be 

\begin{center} \begin{picture} (300 , 100) (0 , 0) 
\GlueArc (150 , 50) (40 , , 180) {5}{8} 
\GlueArc (150 , 50) (40 , 180 , 360) {5}{8} 
\Gluon(50,50) (110,50){5}{4} \Vertex(110,50){2} 
\Gluon(190,50) (250,50){5}{4} \Vertex(190,50){2> 
\end{picture} \\ {\sl A gluon loop diagram} \end{center} 

This code would result in: 




A gluon loop diagram 



The syntax and the meaning of these command are explained in the next section. 
One should note that all coordinates are presented in units of 1 point. There 
are 72 points in an inch. It is possible to use scale transformations if these units 
are not convenient. 
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Currently the primitives are mainly useful for the drawing of Feynman di- 
agrams and the drawing of flowcharts. This means that the commands were 
designed to draw a number of these graphs. Of course many more things can 
be drawn with them, like scatter plots, histograms etc. With such use it may 
be though that the user will miss some primitives that would come in handy. 
Of course in future versions there will be more primitives, especially when users 
suggest some really good ones. 

2 The commands 

The commands that are currently available in Axodraw are (in alphabetic or- 
der): 

• \ Arrow Arc(x,y)(r,0i, 02 ) 

Draws an arc segment centered around (x,y). The radius is r. The arc- 
segment runs counterclockwise from 0i to 02- All angles are given in 
degrees. In the middle of the segment there will be an arrow. 

• \ArrowArcn(x,y) (radius, 0i, 

Draws an arc segment centered around (x,y). The radius is r. The arc- 
segment runs clockwise from 0i to 02- All angles are given in degrees. In 
the middle of the segment there will be an arrow. 

• \ArrowLine(a;i,yi)(a;2,y2) 

Draws a line from (xi,yi) to (£2 ,2/2)- There will be an arrow in the middle 
of the line. 

• \BBox(a;i,yi)(a;2,y2) 

Draws a box of which the contents are blanked out. This means that 
anything that was present at the position of the box will be overwritten. 
The lower left corner of the box is at (x\,y\) and {x2,y2) is the upper right 
corner of the box. 

• \BBoxc(x,y)(width,height) 

Draws a box of which the contents are blanked out. This means that 
anything that was present at the position of the box will be overwritten. 
The center of the box is at (x,y). Width and height refer to the full width 
and the full height of the box. 

. \BCirc(x,y){r} 

Draws a circle of which the contents are blanked out. This means that 
anything that was present at the position of the circle will be overwritten. 
The center of the circle is at (x,y). r is its radius. 

• \Boxc(x,y)(width,height) 

Draws a box. The center of the box is at (x,y). Width and height refer 
to the full width and the full height of the box. 
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• \BText(x,y){text} 

Draws a box with one line of centered postscript text in it. The box is 
just big enough to fit around the text. The coordinates refer to the center 
of the box. The box is like a BBox in that it blanks out whatever was at 
the position of the box. 

• \B2Text(x,y){textl}{text2} 

Draws a box with two lines of centered postscript text in it. The box is 
just big enough to fit around the text. The coordinates refer to the center 
of the box. The box is like a BBox in that it blanks out whatever was at 
the position of the box. 

• \CArc(x,y) (radius, 4>\,4>2) 

Draws an arc segment centered around (x,y). The radius is r. The arc- 
segment runs counterclockwise from (f>\ to (f>2- All angles are given in 
degrees. 

• \Curve{(a;i,yi)(a;2,?/2) ■ ■ ■ (z„,y„)} 

Draws a curve through the given points. The x-values are supposed to 
be in ascending order. The curve is a combination of quadratic and third 
order segments and is continuous in its first and second derivatives. 

• \DashArrowArc(x,y) (r,0i ,</>2){dashsize} 

Draws a dashed arc segment centered around (x,y). The radius is r. The 
arc-segment runs counterclockwise from (f>\ to (f>2- All angles are given in 
degrees. In the middle of the segment there will be an arrow. The size of 
the dashes is approximately equal to 'dashsize'. 

• \DashArrowArcn(x,y) (radius, (^1,(^2) {dashsize} 

Draws a dashed arc segment centered around (x,y). The radius is r. The 
arc-segment runs clockwise from <fii to <f>2- All angles are given in degrees. 
In the middle of the segment there will be an arrow. The size of the 
dashes is approximately equal to 'dashsize'. 

• \DashArrowLine(a;i ,2/1 ) (X2 ,2/2 ) {dashsize} 

Draws a line from (x\,yi) to (£2,2/2) with a dashed pattern. The size of the 
black parts of the pattern is given by 'dashsize'. The alternating pieces 
have equal length. The size of the pattern is adjusted so that both the 
begin and the end are black. Halfway the line there is an arrow. 

• \DashCArc(x,y) (radius, (^1,(^2) {dashsize} 

Draws a dashed arc segment centered around (x,y). The radius is r. The 
arc-segment runs counterclockwise from <fi\ to <f>2- All angles are given in 
degrees. The size of the dashes is determined by 'dashsize'. This size is 
adjusted somewhat to make the result look nice. 

• \DashCurve{(a;i, yi)(x2, 2/2) • • • {x n , y n )}{ dashsize} 

Draws a dashed curve through the given points. The x-values are sup- 
posed to be in ascending order. The curve is a combination of quadratic 
and third order segments. The size of the black parts and the white parts 
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will be approximately 'dashsize' each. Some adjustment takes place to 
make the pattern come out right at the endpoints. 

• \DashLine(2:i, 2/1) (22 ,2/2) {dashsize} 

Draws a line from (2:1,2/1) to (£2,2/2) with a dashed pattern. The size of the 
black parts of the pattern is given by 'dashsize'. The alternating pieces 
have equal length. The size of the pattern is adjusted so that both the 
begin and the end are black. 

• \EBox(a;i,yi) (2:2,3/2) 

Draws a box. The lower left corner of the box is at (xi,yi) and (2:2,2/2) is 
the upper right corner of the box. 

• \GBox(2:i,?/i)(2!2,2/2){grayscale} 

Draws a box. The lower left corner of the box is at (2:1,2/1) and (2:2,2/2) 
is the upper right corner of the box. The contents of the box are lost. 
They are overwritten with a color gray that is indicated by the parameter 
'grayscale'. This parameter can have values ranging from (black) to 1 
(white). 

• \GBoxc(x,y) (width, height) {grayscale} 

Draws a box. The center of the box is at (x,y). Width and height refer 
to the full width and the full height of the box. The contents of the box 
are lost. They are overwritten with a color gray that is indicated by the 
parameter 'grayscale'. This parameter can have values ranging from 
(black) to 1 (white). 

• \GCirc(x,y){radius}{grayscale} 

Draws a circle around (x,y) with radius r. The contents of the circle are 
lost. They are overwritten with a color gray that is indicated by the 
parameter 'grayscale'. This parameter can have values ranging from 
(black) to 1 (white). 

• \GlueArc(x,y)(r,<^i,<^2){amplitude}{windings} 

Draws a gluon on an arc-segment. The center of the arc is (x,y) and r is its 
radius. The arc segment runs counterclockwise from 4>i to fa. The width 
of the gluon is twice 'amplitude', and the number of windings is given by 
the last parameter. Note that whether the curls are inside or outside can 
be influenced with the sign of the amplitude. When it is positive the curls 
are on the inside. 

• \Gluon(2;i,2/i)(2:2,2/2){amplitude}{windings} 

Draws a gluon from (21,2/1) to (2:2,2/2)- The width of the gluon will be 
twice the value of 'amplitude'. The number of windings is given by the 
last parameter. If this parameter is not an integer it will be rounded to 
an integer value. The side at which the windings lie is determined by the 
order of the two coordinates. Also a negative amplitude can change this 
side. 
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• \GOval(x,y)(h,w)(0) {grayscale} Draws an oval with an internal color in- 
dicated by grayscale. This parameter can have values ranging from 
(black) to 1 (white). The center of the oval is given by (x,y). Its height is 
h, and the width is w. In addition the oval can be rotated counterclock- 
wise over 4> degrees. The oval overwrites anything that used to be in its 
position. 

• \GText(x,y) {grayscale} {text} 

Draws a gray box with one line of centered postscript text in it. The box 
is just big enough to fit around the text. The coordinates refer to the 
center of the box. The box is like a BBox in that it blanks out whatever 
was at the position of the box. 

• \G2Text(x,y){grayscale}{textl}{text2} 

Draws a gray box with two lines of centered postscript text in it. The 
box is just big enough to fit around the text. The coordinates refer to the 
center of the box. The box is like a BBox in that it blanks out whatever 
was at the position of the box. 

• \LinAxis(a;i,2/i)(a;2,y2)(-^D Ahashsize , offset, width) 

This draws a line to be used as an axis in a graph. Along the axis are hash 
marks. Going from the first coordinate to the second, the hash marks are 
on the left side if 'hashsize', which is the size of the hash marks, is positive 
and on the right side if it is negative. Njj is the number of 'decades', 
indicated by fat hash marks, and d is the number of subdivisions inside 
each decade. The offset parameter tells to which subdivision the first 
coordinate corresponds. When it is zero, this coordinate corresponds to a 
fat mark of a decade. Because axes have their own width, this is indicated 
with the last parameter. 

• \Line(a;i, 2/1X2:2,2/2) 

Draws a line from (2:1,2/1) to (2:2,2/2)- 

• \LogAxis(2;i, 2/1) (2:2, 2/2) (-^l, hashsize ,offset,width) 

This draws a line to be used as an axis in a graph. Along the axis are 
hash marks. Going from the first coordinate to the second, the hash marks 
are on the left side if 'hashsize', which is the size of the hash marks, is 
positive and on the right side if it is negative. Nl is the number of orders 
of magnitude, indicated by fat hash marks. The offset parameter tells 
to which integer subdivision the first coordinate corresponds. When it 
is zero, this coordinate corresponds to a fat mark, which is identical to 
when the value would have been 1. Because axes have their own width, 
this is indicated with the last parameter. 

• \LongArrow(a;i, 2/1) (2:2,2/2) 

Draws a line from (2:1,2/1) to (2:2,2/2)- There will be an arrow at the end 
of the line. 



6 



• \LongArrowArc(x,y)(r,0i,02) 

Draws an arc segment centered around (x,y). The radius is r. The arc- 
segment runs counterclockwise from fa to fa. All angles are given in 
degrees. At the end of the segment there will be an arrow. 

• \LongArrowArcn(x,y) (radius, (^1,(^2) 

Draws an arc segment centered around (x,y). The radius is r. The arc- 
segment runs clockwise from fa to fa. All angles are given in degrees. At 
the end of the segment there will be an arrow. 

• \Oval(x,y)(h,w)(</>) Draws an oval. The center of the oval is given by 
(x,y). Its height is h, and the width is w. In addition the oval can be 
rotated counterclockwise over cj) degrees. The oval does not overwrite its 
contents. 

• \Photon(rri ,yi ) (2; 2 ,2/2) {amplitude} {wiggles} 

Draws a photon from (xi,y\) to (X2 ,2/2)- The width of the photon will be 
twice the value of 'amplitude'. The number of wiggles is given by the last 
parameter. If twice this parameter is not an integer it will be rounded 
to an integer value. Whether the first wiggle starts up or down can be 
influenced with the sign of the amplitude. 

• \PhotonArc(x,y)(r,<^i,<^2){amplitude}{wiggles} 

Draws a photon on an arc-segment. The center of the arc is (x,y) and r 
is its radius. The arc segment runs counterclockwise from fa to fa. The 
width of the photon is twice 'amplitude', and the number of wiggles is 
given by the last parameter. Note that the sign of the amplitude influences 
whether the photon starts going outside (positive) or starts going inside 
(negative). If one likes the photon to reach both endpoints from the 
outside the number of wiggles should be an integer plus 0.5. 

• \PText(x,y)((^)[mode]{text} 

Places a postscript text. The focal point is (x,y). The text is the last 
parameter. The mode parameter tells how the text should be positioned 
with respect to the focal point. If this parameter is omitted the center 
of the text will correspond to the focal point. Other options are: 1 for 
having the left side correspond to the focal point, r for having the right 
side correspond to it, t for having the top at the focal point and b for the 
bottom. One may combine two letters as in [bl], as long as it makes sense. 
The parameter ^ is a rotation angle. The text is written in the current 
postscript font. This font can be set with the SetPFont command. 

• \rText(x,y) [mode] [rotation] {text} 

Places a rotated text. The focal point is (x,y). The text is the last 
parameter. If the rotation parameter is the character 1 the text will be 
rotated left by 90 degrees, if it is an r it will be rotated to the right by 
90 degrees and when it is the character u the text will be rotated by 
180 degrees. When there is no character there is no rotation and the 
command is identical to the Text command. The mode parameter tells 
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how the resulting box should be positioned with respect to the focal point. 
If this parameter is omitted the center of the box will correspond to the 
focal point. Other options are: 1 for having the left side correspond to 
the focal point, r for having the right side correspond to it, t for having 
the top at the focal point and b for the bottom. One may combine two 
letters as in [bl], as long as it makes sense. 

• \SetPFont{fontname}{fontsize} 

Sets the postscript font to a given type and scale. 

• \SetScale{scalevalue} 

Changes the scale of all graphics operations. Unfortunately it does not 
change the scale of the text operations (yet?). A 'scalevalue' of 1 is the 
default. It is allowed to use floating point values. 

• \SetOffset(x_offset,y_offset) 

Adds the offset values to all coordinates at the T^X level. This makes it 
easier to move figures around. 

• \ Set ScaledOffset (x_offset ,y .offset ) 

Adds the offset values to all coordinates at the postscript level. This is 
done after scaling has been applied. Hence one can work with the scaled 
coordinates. This can be very handy when drawing curves. 

• \SetWidth{widthvalue} 

Changes the linewidth in all graphics operations. It does not change the 
linewidth of the text operations. That is a matter of font selection. A 
'widthvalue' of 0.5 is the default. It is allowed to use floating point values. 

• \Text(x,y)[mode]{text} 

Places a text. The focal point is (x,y). The text is the last parameter. 
The mode parameter tells how the text should be positioned with respect 
to the focal point. If this parameter is omitted the center of the text will 
correspond to the focal point. Other options are: 1 for having the left side 
correspond to the focal point, r for having the right side correspond to it, 
t for having the top at the focal point and b for the bottom. One may 
combine two letters as in [bl], as long as it makes sense. 

• \Vertex(x,y){r} 

Draws a fat dot at (x,y). The radius of the dot is given by r. 

• \ZigZag(rri,yi)(a;2,y2){amplitude}{wiggles} 

Draws a zigzag line from (x\,yi) to (2:2,2/2)- The width of the zigzagging 
will be twice the value of 'amplitude'. The number of zigzags is given by 
the last parameter. If twice this parameter is not an integer it will be 
rounded to an integer value. Whether the first zigzag starts up or down 
can be influenced with the sign of the amplitude. 
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3 Examples 



Although the previous section contains all the commands and their proper syn- 
tax a few examples may be helpful. 

3.1 Text modes 

The meaning of the mode characters in the text command can best be demon- 
strated. The statements 

\begin{center} \begin{picture} (300 , 100) (0,0) 
\CArc (50 , 75) (2 , , 360) \Text (50 , 75) [It] {left-top} 
\CArc (50 , 50) (2 , , 360) \Text (50 , 50) [1] {left-center} 
\CArc (50 , 25) (2 , , 360) \Text (50 , 25) [lb] {left-bottom} 
\CArc(150,75) (2,0,360) \Text(150,75) [t] {center-top} 
\CArc(150,50) (2,0,360) \Text(150,50) [] {center-center} 
\CArc(150,25) (2,0,360) \Text(150,25) [b] {center-bottom} 
\CArc (250 , 75) (2 , , 360) \Text (250 , 75) [rt] {right-top} 
\CArc (250 , 50) (2 , , 360) \Text (250 , 50) [r] {right-center} 
\CArc (250 , 25) (2 , , 360) \Text (250 , 25) [rb] {right -bottom} 
\end{picture} \end{center} 

produce 9 texts and for each the focal point is indicated by a little circle. It 
looks like 

left-top center-top right-top° 

4eft-center centerQcenter right-centero 
deft-bottom center-bottom right-bottomo 



This illustrates exactly all the combinations of the mode characters and what 
their effects are. The text can of course contain different fonts, math mode and 
all those little things that are usually easier in I^Tf^X than in postscript. 

3.2 The windings of a gluon 

Gluons are traditionally represented by a two dimensional projection of a he- 
lix. Actually close inspection of some pretty gluons reveals that it is usually 
not quite a helix. Hence the gluons in Axodraw are also not quite helices. In 
addition one may notice that the begin and end points deviate slightly from 
the regular windings. This makes it more in agreement with hand drawn glu- 
ons. When a gluon is drawn, one needs not only its begin and end points but 
there is an amplitude connected to this almost helix, and in addition there are 
windings. The number of windings is the number of curls that the gluon will 
have. Different people may prefer different densities of curls. This can effect 
the appearance considerably: 
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\begin{center} 

\begin{picture} (330 , 100) (0 , 0) 
\Gluon(25,15) (25,95){5}{4} 
\Gluon(95,15) (95,95){5}{5} 
\Gluon(165,15) (165,95){5}-[6} 
\Gluon(235,15) (235,95){5}{7} 
\Gluon(305,15) (305,95){5}{8} 
\end{picture} 
\end{center} 



\Text(25,7) [] {4 windings} 
\Text(95,7) [] {5 windings} 
\Text (165,7) [] {6 windings} 
\Text (235,7) [] {7 windings} 
\Text (305,7) [] {8 windings} 



This code results in: 



4 windings 5 windings 6 windings 7 windings 8 windings 



The influence of the amplitude is also rather great. The user should experiment 
with it. There is however an aspect to the amplitude that should be discussed. 
For a straight gluon the amplitude can determine on which side the curls are. 
So does the direction of the gluon: 

\begin{center} 

\begin{picture} (325 , 100) (0 , 0) 
\Gluon(50,15) (50,95){5}{6} 

\Text(50,7) []{amp $> 0$} \Text (40, 50) [] {$\uparrow$} 
\Gluon(125,95) (125,15){5}{6} 

\Text(125,7) []{amp $> 0$} \Text (115 , 50) [] {$\downarrow$} 
\Gluon(200,15) (200,95){-5}{6} 

\Text(200,7) []{amp $< 0$} \Text (190, 50) [] {$\uparrow$} 
\Gluon(275,95) (275,15){-5}{6} 

\Text(275,7) []{amp $< 0$} \Text (265 , 50) [] {$\downarrow$} 

\end{picture} 

\end{center} 

The picture gets the following appearance: 




amp > amp > amp < amp < 
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For straight gluons one does not need the option of the negative amplitude. 
It is however necessary for gluons on an arc segment. In that case the arc is 
always drawn in an anticlockwise direction. Hence the direction is fixed and 
only the amplitude is left as a tool for determining the side with the curls. 

3.3 Scaling 

Sometimes it is much easier to design a figure on a larger scale than it is needed 
in the eventual printing. In that case one can use a scale factor, either during 
the design or in the final result. We use the figure in the first section as an 
example: 

\vspace{-10pt> \hfill \\ 
\SetScale{0.3} 
\begin{picture} (70 , 30) (0 , 13) 
\GlueArc (120 , 50) (40 , , 180) {5}{8} 
\GlueArc (120 , 50) (40 , 180 , 360) {5}{8} 
\Gluon(20,50) (80,50){5}{4> \Vertex(80,50){2> 
\Gluon(160,50) (220,50){5}{4} \Vertex(160,50){2} 
\end{picture} $+$ others 

$ = C_A(\frac{5H3}+\frac{3lH9}\epsilon) 
+ n_F(-\frac{2}{3}-\frac{10}{9}\epsilon)$ 
\vspace{10pt> \hfill \\ 

We have lowered the figure by 13 points (the (0,13) in the picture statement) 
to make it look nice with respect to the equal sign. The result is 



This way it is rather straightforward to make whole pictorial equations. Of 
course some things are not scale invariant. The appreciation of a figure may 
be somewhat different when the scale is changed. In the above case one might 
consider changing the amplitude of the gluons a little bit. Changing this from 
5 to 7 and at the same time reducing the number of windings from 4 to 3 for 
the straight gluons and from 8 to 7 for the gluons in the arcs gives 



At this scale this may please the eye more. 

There is one problem with scaling. Currently it is only possible to have text 
scale with the rest of a figure when the text has been printed with the PText 
command. This makes the typesetting more complicated, but the scaling of the 
TgX pixel fonts would give rather poor results anyway. 

3.4 Photons 

When drawing photons one should take care that the number of wiggles is 
selected properly. Very often this number should be an integer plus 0.5. This 
can be seen in the following example: 




+ others = C A {\ + f e ) + «f (~| " T e ) 




+ others = CU(§ + f-e) + n F (-| - ^e) 
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\begin{center}\begin{picture} (300 , 56) (0 , 0) 
\ArrowLine(100,10) (200,10) 

\PhotonArc(150,10) (30,0,180){4}{8.5} °/„ 8.5 wiggles 
\LongArrowArc ( 150 , 10) (20 , 60 , 120) 
\Vertex(180,10){1.5} \Vertex(120 , 10) {1 . 5} 
\end{picture} \end{center} 

This gives the 'proper' picture as it would usually drawn by hand: 




When the number of wiggles is reduced to 8 we obtain: 




This is not as nice. Somehow the symmetry is violated. One should also take 
care that the wiggles start in the proper way. If we make the amplitude negative 
we see that the photons are not 'right' either: 




Sometimes these things require some experimenting. 
3.5 Flowcharts 

There are several commands for creating boxes with text in them. This can be 
a box with either one line of text or with two lines of text. The rest is just a 
matter of drawing lines and circle segments with arrows. If the text is to scale 
with the picture one needs to use the postscript fonts. The result of scaling the 
TgX fonts is usually rather ugly, because these fonts are pixel fonts. Here we 
present an example. It might describe a system for the automatic computation 
of cross-sections: 

\begin{center} \begin{picture} (320 , 320) (0,0) 

\SetPFont{Helvetica}{10} 

\SetScale{0.8} 

\ArrowLine (200 , 40) (200 , 10) \ArrowLine (200 , 100) (200 , 40) 
\ArrowLine(200,150) (200,100) \ArrowLine(100, 130) (200,100) 
\ArrowLine (85 , 95) (200 , 100) \ArrowLine (260 , 105) (200 , 100) 
\ArrowLine(250,135) (200,100) \ArrowLine(160,75) (200,100) 



12 



\ArrowLine (200 , 100) (250 , 70) \ArrowLine (200 , 185) (200 , 150) 
\ArrowLine (200 , 220) (200 , 185) \ArrowLine (200 , 250) (200 , 220) 
\ArrowLine (240 , 263) (200 , 250) \ArrowLine (240 , 237) (200 , 250) 
\ArrowLine (200 , 285) (200 , 250) \ArrowLine (200 ,310) (200 , 285) 
\ArrowLine(200,335) (200,310) \ArrowLine(180,360) (200,335) 
\ArrowLine(200,385) (180,360) \ArrowLine(50,370) (180,360) 
\ArrowArc (200 , 247 . 5) (62 . 5 , 90 , 180) 
\ArrowArc (200 , 247 . 5) (62 . 5 , 180 , 270) 

\ArrowLine(210,385) (300,360) \ArrowLine(210,335) (300,360) 
\ArrowLine(80,300) (80,130) \ArrowLine(190,335) (80,300) 
\ArrowLine(190,385) (80,300) \ArrowLine(50,335) (80,300) 
\ArrowLine(300,360) (340,340) \ArrowArcn(205,347.5) (37.5,90,270) 

\BCirc (200 , 100) {10} \BCirc (200 , 100) {5} 

\BCirc (200 , 40) {7 . 5} \BCirc (200 , 250) {10} 

\BCirc (200 , 250) {5} \BCirc (200 ,310) {7 . 5} 

\BCirc (180 , 360) {7 . 5} \BCirc (80 , 300) {7 . 5} 

\BCirc (300 , 360) {7 . 5} \GCirc (200 , 185) {7 . 5}{0 . 9} 

\BText (200, 285) {Form program} \BText (200, 335) {Diagrams} 

\BText (200 , 385) {Model} \BText (200 , 10) {events} 

\BText (80 , 95) {Axolib} \BText (350 , 335) {Pictures} 

\GText (137 . 5 , 247 . 5) {0 . 9}{instructions} 

\B2Text (260 , 70) {Cross-sections}{Histograms} 

\B2Text (140 , 75) {Monte Carlo}{Routine} 

\B2Text(275,105){FF}{l loop integrals} 

\G2Text (260 , 135) {0 . 9}{Spiderlib}{Fortran/C} 

\G2Text (200 , 150) {0 . 9}{Matrix}{Element} 

\B2Text (80 , 130) {Kinematics}{Conf iguration} 

\G2Text (200 , 220) {0 . 9}{0utput}{Formula} 

\G2Text (260 , 263) {0 . 9}{Spiderlib}{Form part} 

\G2Text (260 , 237) {0 . 9}{FF supportHlibrary} 

\B2Text (40 , 370) {Reaction}{selection} 

\B2Text (40 , 340) {Specif ication}{Cuts , etc . } 

\PText(211,36) (0) [lb] {Event Generator} 
\PText (211 , 181) (0) [lb] {Code Generator} 
\PText ( 162 , 258) (0) [lb] {FORM} 

\PText(211,301) (0) [lb] {Form program construction} 
\PText(191,362) (0) [lb] {Diagram} 
\PText(191,352) (0) [lb] {Generator} 
\PText (31 1,370) (0) [lb] {Postscript} 
\PText (31 1,360) (0) [lb] {Generator} 
\PText (91 , 292) (0) [lb] {Kinematics} 
\PText (91 , 282) (0) [lb] {Generator} 
\end{picture} \end{center} 

This gives the chart 
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(j) Ev ent Generator 
events 



3.6 Curves and graphs 

Axodraw is equipped with a curve fitting facility that can draw smooth curves 
through a set of coordinates. Coupled to this is a set of commands to draw 
the axes that are typically needed for the use of graphs and histograms. An 
example of a complete picture would be 

\begin{center} \begin{picture} (360 , 440) (0,0) 
\Set0ff set (40,30) 

\LinAxis(0,0) (300,0) (3,10,5,0,1.5) 

\LinAxis (0 , 400) (300 , 400) (3 , 10 , -5 , , 1 . 5) 

\LogAxis(0,0) (0,400) (4,-5,2,1.5) 

\LogAxis (300 , 0) (300 , 400) (4,5,2,1.5) 

\SetScale{100.} \SetWidth{0 . 005} 

\Curve{ ( . 1057001 , 1 . 2997) ( . 1057003 , 1 . 5399) 

( . 1057006 , 1 . 6908) ( . 1057010 , 1 . 8019) ( . 1057030 , 2 . 0406) 

( . 1057060 , 2 . 191 1) ( . 1057100 , 2 . 3020) ( . 1057300 , 2 . 5403) 

( . 1057600 , 2 . 6904) ( . 1058000 , 2 . 8007) ( . 1060000 , 3 . 0365) 

( . 1080000 , 3 . 4512) ( . 1 100000 , 3 . 5600) ( . 1200000 , 3 . 6950) 

( . 1300000 , 3 . 6969) ( . 1500000 , 3 . 6308) ( . 1800000 , 3 . 5024) 

( . 2200000 , 3 . 3413) ( . 3000000 , 3 . 0788) ( . 5000000 , 2 . 6374) 

(.8000000,2.2295) (1.0000000,2.0357) (1.3000000 

,1.8078) (1.6000000,1.6275) (2.0000000,1.4336) 
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(2 . 5000000 , 1 . 2398) (3 . 0000000 , 1 . 0815) } 

\DashCurve{( 1.7853600, .0111) (1.7853800, .0228) 

(1.7854000, .0339) (1 . 7856000 ,. 1218) (1.7860000, .2324) 

(1.7870000, .3821) (1 . 7900000 ,. 5786) (1.8000000, .8089) 

(1.8200000, .9765) (1.8500000,1.0869) (1.9000000 

, 1 . 1718) (2 . 0000000 , 1 . 2335) (2 . 1000000 , 1 . 2468) 

(2 . 2000000 , 1 . 2413) (2 . 4000000 , 1 . 2064) 

(2 . 7000000 , 1 . 1340) (3 . 0000000 , 1 . 0574) }{0 . 05} 

\SetScale{l . }\SetWidth{0 . 5} 

\Line(200,360) (270,360) \Text(195,360) [r]{ 

\large$e~+e~-\rightarrow\mu~+\mu~-$} 

\DashLine(200,330) (270,330){5}\Text(195,330) [r]{ 

\large$e~+e~-\rightarrow\tau~+\tau~-$} 

\Text(0,-10) []{0} \Text(100,-10) []{1} 

\Text(200,-10) []{2> \Text(300,-10) []{3> 

\Text(150,-25) [] {\large Beam energy in GeV} 

\Text (-10,70) []{$1$> \Text(-10,170) []{$10$} 

\Text(-10,270) []{$10"2$} \Text(-10,370) []{$10"3$} 

\rText(-25,220) [] [1] {\Large$\sigma$ in nb} 

\ArrowLine(190,270) (160,300) 

\ArrowLine(160,240) (190,270) 

\ArrowLine (270 , 300) (240 , 270) 

\ArrowLine (240 , 270) (270 , 240) 

\Photon( 190 , 270) (240 , 270) {4}{4 . 5} 

Wertex (190 , 270) {1 . 5} \ Vertex (240 , 270) {1 . 5} 

\end{picture} \\ {\sl \hskip 10 pt Threshold 

effects for $\mu$ and $\tau$} \end-fcenter} 

and the resulting picture would be 
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12 3 

Beam energy in GeV 
Threshold effects for [x and r 

Of course one can scale these pictures further, but because the scale factor has 
been used to enter the data points these should then be adapted too. Note that 
when the scale is blown up by a factor 100, the linewidth has to be scaled down 
or disasters will take place. 
Finally a playful example: 

\begin{center}\begin{picture} (300 , 56) (0 , 0) 
\Gluon(150,25) (200,25){3}{6} 
\Photon(150,35) (200,45){3}{6} 
\ZigZag(150,15) (200,5){3}{6} 
\Line(100,25) (150,25) 
\G0val(150,25) (20,10) (0){0.5> 
\end{picture} \end{center} 
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which results in 



Acknowledgement: The author wishes to thank G.J. van Oldenborgh for help 
with some of the T^X macros. 

Axodraw can be obtained by means of anonymous ftp from ftp.nikhef.nl. It 
is located in the directory pub/form/axodraw. Commentary and suggestions 
should be sent to the author at t68@nikhef.nl. 
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